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, Abstract 

m ■ 

^ ' We generalize Israel's formalism to cover singular shells embedded 

' in a non- vacuum Universe. That is, we deduce the relativistic equa- 

. tion of motion for a thin shell embedded in a Schwarzschild/Friedmann- 

' Lemaitre-Robertson-Walker spacetime. Also, we review the embedding of 
a Schwarszchild mass into a cosmological model using "curvature" coor- 

' dinates and give solutions with (Sch/FLRW) and without the embedded 

g ■ mass (FLRW). 

I . 1 Introduction 

^ . 

• • , The first applications of Israel's relativistic theory of singular surfaces |P |2] 

, ^ ' were concerned with neutral and charged surfaces layers in static spacetimes 

^ , with or without a cosmological constant [H] - 0. Surface layers with different 

^ ' equations of state, i.e. different relationships between pressure and density were 

5^ , considered, for example domain walls Also, methods to quantize a shell in 

general relativity has been developed and studied for static spacetimes pi I10|. 
Lately Israel's theory has been applied to three-dimensional branes in a five- 
dimensional bulk [111 I12j . 

In the present work we want to extend the application of Israel's formalism 
to cases with singular layers (two-dimensional branes) in expanding Universe 
models. In the static case with a neutral surface layer there is Schwarzschild 
or Schwarzschild-de Sitter spacetime. For non-static spacetimes the relativistic 
equation of motion for surface layers has been given for Friedmann-Lemaitre- 
Robertson- Walker (FLRW) spacetimes ^S] - QEj- In the more general case we 
must consider spacetime outside a mass embedded in an expanding Universe 
model. A description of this was developed by R. Gautreau [T7I[TH| . 
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We shall therefore apply Israel's theory to a singular layer with the Gautreau 
metric outside the shell. The equation of motion for the general case is deduced 
here. 

A similar problem, with a slowly rotating shell in an almost FLRW dust 
universe model was considerd by Klein |19j without using Israel's formalism 
and later by T. Dolezel, J. Bicak and N. Deruelle [201 using Israel's formalism. 
They restricted their treatments by considering only a shell comoving with the 
cosmic dust outside the shell, and with Minkowski spacetime inside the shell. 
In our treatment we consider a non-rotating shell that needs not be comoving 
with the cosmic fluid. 



2 Schwarzschild mass embedded into a cosmo- 
logical model 

Embedding a Schwarzschild mass, m, into a cosmological model is most easily 
done in "curvature" coordinates. That is, in coordinates for which the radial 
coordinate R gives an angular part R^dil^. The embedding of a Schwarzschild 
mass into a spatially flat FLRW Universe model was given by R. Gautreau 
in Also, a thorough investigation of the FLRW models with vanishing 

cosmological constant, A = 0, in these coordinates has been performed by R. 
Gautreau, see |21j and references therein. In this section we give a short review. 
Then we give the equations for embedding a Schwarzschild mass into a general 
FLRW Universe model and provide solutions with and without the embedded 
mass. 

2.1 The metric 

For a spherically symmetric spacetime we can write a general metric as: 

ds' = T)f{R, T)dT' + + R'dn', (1) 

where 

dn^ ^ de^ + siu^ ed<j)'^ . (2) 

Here A and / are functions of T and R that are settled by Einstein's field equa- 
tions. The physical interpretation of the time coordinate T is that it measures 
the time on clocks that are located at points for which R = constant relative 
to our chosen origin R = 0} In the FLRW Universe models one can define a 
global time coordinate. Hence, we wish to record time on a geodesically moving 
clock. To transform to a geodesically moving clock we consider the geodesic 
equation, V^j^V^ = 0, for a radially moving clock, i.e d6 = d(j) = 0; where is 
the tangent vector to the geodesic curve. For timelike geodesies, V^ Vf^ = —1,^ 
the solution is found to be (see |17|): 

^The time laps in T and the time laps recorded on clocks that measures proper time t 
(ds^ = —dr'^) for i?,, 9,<f> = constant are related by dT = ^^^ dr. 

^We use units where the speed of light and the gravitational constant are set to unity, 
c = G = 1. 
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V'^{T,R)^^^(^-^,,{£'-A)K0,0^ , , = ±1, (3) 
with the conditions (coordinate transformations): 

t,n = -U-\£'-A)i (4) 

t,T = £f, (5) 



where £ is an energy parameter for the reference particles, i.e. the geodesic 
clocks. £ depends on our choice of reference system, and can be used to de- 
scribe the open, closed and flat FLRW Universe models. The time coordinate 
t measures the time recorded on clocks moving on radial geodesies, i.e. t mea- 
sures proper time, t = t, and thus on radial geodesies we have ds^ — —dt^. The 
sign given by <j indicates whether the geodesic clock is moving with increasing 
R (iT = +1) or with decreasing R (<r = — 1). So that for an expanding Universe 
we have = Now we can make a coordinate transformation from (T, R) to 
(t, R) coordinates. In (t, R) coordinates the radial geodesies are described by 
the 4-velocity 

V^{t,R) = (1, 0, 0). (6) 

The resulting form of the line element is: 

ds^ = £-^(^- Adi^ + - 2^{£^ - A)^dRdt^ + R^dn"^ 

= --dt^ + £-'^(dR-^{£'^ - Af^dt^ + R^dn^. (7) 

Because t measures time on clocks moving relative to i? = constant the metric 
is non-diagonal. The cosmic particles, e.g. the galaxies are assumed to follow 
radial geodesies. Without the embedded Schwarzschild mass, i.e. m = 0, we 
require the Universe to be isotropic and homogeneous. The choice of coordinates 
ensures the isotropic condition. For a flat Universe model the homogeneity 
condition settles £ ^7]. First we note that for m — Q the coordinates are 
required to reduce to the Minkowskian form at i? = 0, i.e. A = f = 1. Then 
from ^ we find that at i? = the energy parameter £ is given by the Lorentz 
factor 7: f = [l — v^^)^^, where is the coordinate velocity, of the reference 
particles at i? = 0; and because the space-time metric here is flat this is the 
velocity measured by an observer at i? = 0. Hence, £ gives the energy to rest 
mass ratio for the cosmic particles. Particles with £ <1 will never reach R — Q, 
thus making a hole in the cosmic fluid. While for £ > 1 they will have some 
velocity, > 0, at i? = 0. This makes i? = a source (<; = -1-1) or a sink 
(<r = — 1) for cosmic particles. For a homogeneous Universe i? = cannot have 
any special significance. Hence, we are lead to 

£^l (8) 

for a spatially flat Universe. If we embed a Schwarzschild mass into this model 
we get an inhomogeneous model, but we still have £ = 1 !17l. Gautreau argues 
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that for a vacuum Universe, i.e. the de Sitter model to have a physical reference 
system we cannot set £ ^ I. The vacuum Universe in the non-diagonal metric 
(|7|) is discussed in For a thorough investigation of the kinematics of the 

de-Sitter Universe see 

For a flat Universe model, i.e. £ = 1, the surface of simultaneity dt = in 
{T)) gives the flat metric. 

dl^ ^ dR^ + R^dfl^. (9) 

In this case R measures the actual distance between the cosmic particles. 

We shall now consider £ for open and closed models. To get a feeling for the 
coordinates in (0 we relate them to the commonly used comoving coordinates 
in cosmology. In these coordinates the metric has the form 

ds^ = -dt^ + a{tf ( + r^nA (10) 
\ 1 — kr-^ J 

where a(t) is the expansion factor and the radial coordinate r is constant along 
the trajectory of the cosmic particles, e.g. galaxies, and t is the time measured on 
clocks moving with them. The sign k gives the spatial topology (i.e. the global 
structure of the t = constant surface): fc = -1-1 leads to a spherical (closed) 
geometry, k ~ —\ gives a hyperbolic (open) geometry and k = Q describes an 
Euclidean (flat) geometry. 

The time coordinates in ^ and (|10|l measure time on clocks moving with 
the cosmic particles. Thus we identify the time coordinates in |2J) and H10|l . 
The radial coordinates are related hy R — ar and the energy parameter, £, for 
the reference particles is given hy £^ = 1 — kr^ = 1 — kC{t, R), where C{t, R) is 
a function that is constant along the radial geodesies, i.e. £ must be constant 
along the geodesies. The Einstein equations give this function up to a constant 
factor. The Hubble factor, H = is now given by 

f.H«^ ,11, 

2.2 The Einstein equations for curvature coordinates 

The Einstein equations with cosmological constant G^^ + A(5^^ = SttT'^ which 
we will need to deal with are: 



-8ttR^T\ + Ai?2 (12) 
SnR^T^ (13) 



d{R{l- 


A)) 


dR 




d{R{l~ 


A)) 


dt 




d{R{l- 


A)) 


dR 



These equations are valid in coordinates a;^ = (t, R, 9, (/>) for an arbitrary energy- 
momentum tensor T'^. The field equations involving second-order derivatives 
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Cpsiaic fluid 



Figure 1: Schwarzschild mass embedded in a cosmic fluid. 



are the angular part, for which T g — T these equations are contained in the 
divergence of the Einstein tensor, Gf'.^ = 0. 

To proceed an energy-momentum tensor for the Universe with an embedded 
Schwarzschild mass is assumed. One may imagine the galaxies to be particles 
of a fluid and that this cosmic fluid fills the whole spacetime.^ The embedded 
Schwarzschild mass, m, is placed with its center at i? = 0, and with boundary 
at R = Rb, see figure As mentioned above the particles of the cosmic fluid 
are assumed to follow the radial geodesies V^. The cosmic fluid is described as 
an ideal fluid. Thus, outside Ri, the energy-momentum tensor is: 

T^i = ip + p)V^V,+pS^^, R>Ri„ (15) 

where p is the mass density and p the pressure of the cosmic fluid. Inserting 
given in © into the energy-momentum tensor (|15|) gives the components: 



(16) 

T*^ = 0, T^^-,i£'-A)iip + p). 

Inside Rf, we will assume that the tt component of the energy-momentum 
tensor can be written as 

T\ = -P + P\, R<Rb, (17) 

^The embedding of a Schwarzschild mass into a cosmological model has also been considered 
by Einstein and Straus 24 (the Swiss Cheese model). They looked at a cosmic dust fluid with 
a vacuum region (A = 0) with a Schwarzschild mass at the center. And showed that one can 
join the Schwarzschild metric smoothly onto the cosmic metric IIUI . An explicit form of this 
metric was found by Schiicking 25_ using "curvature" coordinates. For recent developments 
on this embedding see e.g. 1191 1261 1271 . 
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where p is the energy density of the cosmic fluid and p*j is the part of the energy- 
momentum tensor giving the energy density of the embedded Schwarzschild 
mass. If we assume that p*^ is time independent this is the only component of 
the energy momentum tensor inside Rb we will need. The constant embedded 
Schwarzschild mass, m, bounded by Rb is defined by: 

rRb 

m = -47r / p\R^dR. (18) 



We note that setting m = leads to a description of cosmology in "curvature" 
coordinates. Let us also define a mass function M{t, R) for the Universe: 

M{t, R)^Att I p{t, R)R^dR, (19) 

where the integration is over a t = constant surface. In terms of M the density 
is given by 

^M^^ = A^R^p{t,R), (20) 

where AttR^ gives the area of a sphere centered on i? = 0. 

In the case where we have an embedded mass m the density p of the cosmic 
fluid will have some R dependence and we cannot bring p outside the integrals. 
But for m = the cosmic fluid is homogeneous, and p is only a function of t 
and thus can be brought outside the integrals. 

We now move on to solve the Einstein equations (|12I14|I outside Rb^ with an 
energy-momentum tensor given by (fTH|l and ifTTI . We start by integrating ifT^ 
over a t = constant surface. Using (|18|) and (|19|l this gives the metric function 
in the form 

2(m + Af) A , 
A=l- ^ ^ -3^' (21) 

with the metric given in 0. Note that in general M is a function of t and R. 
We need to determine M. Inserting (|21|l into (|13|l gives 



i?2 7o dt ' 

which gives the energy (mass) fiux across a R — constant surface. The 
component of the energy momentum tensor for an ideal fluid is given in 116|l . 
Equating these two expressions for and using 12()|l we find 

^+<^i£'-A)i^ = ^,4nRMS'-A)K (23) 

which is the partial differential equation determining M. From l|23(l and H2()|l . 
assuming p > —p, we see that for an expanding Universe {<; ~ we have 
< 0; and for a contracting Universe (<; = — 1) we get > 0. 
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To proceed from here we assume an equation of state p — p{p) for the cosmic 
fluid. We assume that the cosmic fluid obeys the barotropic equation of state 



p^up, w^Lo{t,R), (24) 

e.g. uj — Q gives a dust Universe, uj = ^ describes a radiation dominated 
Universe and lj ^ —\ represents a vacuum Universe. Using H2(J|) we get 

— +^(l+^)(f2-A)^— = 0. (25) 
Writing l|23|) in terms of p jJT] for a general Universe leads to 



i?^^ ^R'^dR + <;R^{p+p)(^£^-l + 8TT pR^dR+2m + ^R^y = 0. (26) 

Evaluating along the geodesies we find 

dM , ^, dR 

For dust, p — this gives = 0. So that M is constant along the geodesies. 

Consider the vacuum case. For Lorentz invariant vacuum energy the energy 
momentum tensor is proportional to the metric |28j . i.e. T^^ = Pj^g^u with the 
vacuum energy density constant. Thus, this type of vacuum energy can be 
incorporated into the cosmological constant A; we have A = Snp^. So that here 
the embedding of a Schwarszchild mass into the de Sitter Universe is given by 
M = 0. The trivial solution of (|25|l is M — constant, and if this is a global 
solution then from the definition of M in (|f 9|l and H2()() this constant must be 
zero. Let now be part of the energy momentum tensor, i.e. included in M. 
For vacuum, uj = —1, the second term in (|25|l or H26(l vanish. This gives = 
showing that M is time independent. For a homogeneous space p^ must then be 
constant. Let us also briefly look at the Einstein equations using the diagonal 
metric . The above equations become for an arbitrary T'^ : 



A{R, T) = 1 + ^ r T^R^dR - ^R^, (28) 
' "^^R^dR, (29) 



i?2 ./n dT 



Inf = 8^^ ^{T\ ~ T^T)dR. (30) 

The integration is now over a T — constant surface. The coordinate transfor- 
mation between the coordinate systems is given by (@J) and From (|30|) we 
see that all vacuum space-times have / = 1. And for a Schwarzschild mass 
embedded in a de-Sitter Universe the metric 1^ is given by A = 1 — ^ — ^R^- 
Let us now consider £. From the Einstein equation (|f 4|l and demanding that 
(|^ is satisfied we find 
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For geodesies this gives ^ = 0, which shows that the Einstein equations require 
that £ is a constant along the streamhnes of the cosmic fluid outside Ri,. For a 
flat Universe model, with and without the embedded Schwarzschild mass, this 
constant must be set to f = 1 . For a Universe with a non zero spatial curvature 
we see from eq. H25|l that £ can be written in terms of M for a dust Universe. 
To recover standard cosmology we require that £^ — 1 is proportional to M. 
Hence, we can write 



2M 
Ri ' 



(32) 



where Ri = constant. This constant depends on which galaxy we wish to follow 
and thus it determines M. E.g. for a closed (fc = +1) Universe Ri represents 
the maximum radius reached by a given galaxy. For further discussion of this 
and the open case see 

Next we will consider the velocity and the acceleration of the cosmic particles. 
From the geodesies in © the coordinate velocity is in general given by 



(33) 



with the Hubble parameter given hy H — Inserting for A and £ gives 

2M 2{m + M) A 



dR 
dt 



Ri 



R 



(34) 



where M is constant for a dust {p — 0) Universe. Differentiating (|34|l and 
inserting (|27|l with p = we get the acceleration for a dust Universe, 



d'^R m + M 

For m = we have a homogeneous cosmic fluid, thus we can find M from 
considering the geodesies H34II and (|35|l . If m 7^ the cosmic fluid is inhomoge- 
neous and we must find M by treating the coordinates, t and R, as independent 
variables, i.e. from eq. (|23|l . Also, we note that for m — equations H34() and 
(|35|l are equivalent to the Friedmann equations, i.e. the dynamical equations 
obtained from Einsten's field equations when using the comoving coordinates in 
For A = (|35|l is identical to Newton's gravitational law, see j29j . 
For light the trajectories are given by ds^ — 0, thus the paths of radially 
moving light are described by 



^^{£'-A)i±£. (36) 

For a closed Universe model, k — +1, in the limit Ri 2M we have £^ ^ 
giving a singularity in the metric (0). In this case 1)34(1 approaches ()36|l . i.e. 
the galaxies approach the speed of light. This is the maximum size of the 
Universe 1291. 
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3 Solutions 



3.1 m = 

For a general Universe without the embedded Schwarzschild mass we can bring 
p outside the integrals above. The mass function, M , for the Universe defined 
in l|19|l becomes 

M = -R\{t). (37) 

For a dust Universe this gives the mass of the Universe inside R at time t. In the 
following we shall consider a flat Universe model, £ — \., such that R measures 
the proper distance between the cosmic particles. Eq. then reduces to: 

f = -^(p + p)(247rp + 3A)i (38) 
Along the geodesies (see (|2Z|)) eq. (p?^ reads 

independent of A (and £). From H38|l and H39|l we find a 1-parameter family 
of solutions. We assume an equation of state of the form p = top with w = 
constant. Integration of then gives p{t) in terms of the geodesies, R{t): 

p[t) = bR{t)-'^'^^+^\ (40) 

This expression for the density is in fact valid for all the FLRW Universe models, 
i.e. regardless of A and S. For dust u) ~ Q which inserted into (|40|l gives 
p (X , i.e. M is constant, for radiation we have ui = ^ and we get p cx 
For vacuum, w = — 1, we get p = constant which is also obtained from eq. (|38|l : 
the expansion of the vacuum Universe is found from eq. H34|l : 



R(xcosh^^t, £<1, (41) 
iiocev^*, £ = 1, (42) 

R(xsmh^t, £>l. (43) 

In the following we assume that w 7^ — 1 . We note that the density p does not 
depend on the integration constant b. This constant depends on which geodesic 
(galaxy) we wish to follow, (see below). Thus, we can relate b to the comoving 
coordinate r in (|10|l . Normalizing the expansion factor, a, so that ao — 1, where 
the index refers to the present time, we get b = r^^^~^^^ po- For a dust Universe 
M is constant and we find that M = ^b. That is, 6 is a measure of the mass 
inside a comoving radius r = constant. Alternatively, one can regard 6 as a 
measure of the energy of a cosmic particle relative to i? = 0, i = 0. 

Solving eq. H38I) for A = gives 

(1 + cj)267rt2 ■ ^"^^^ 
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(a) Einstein-de Sitter 



(b) Friedmann-Lemaitre 



Figure 2: Graphs showing the radial timehke geodesies for different values of b 
in a Universe containing (a) dust, and (b) dust and vacuum energy. 



We have set the integration constant to zero so that the big bang is placed 
at R{t = 0) = 0. Eq. H44|l shows that the density of all ideal fluids, in a flat 
Universe, with p — ujp are proportional to t^^ . Combining (|40|l and (|44|l we 
find the geodesies: 



2\ 3(1+1^) 



(45) 



The geodesies for for the Einstein-de Sitter model, a; = 0, are displayed in fig|21 
where the geodesic i? = represents the trajectory of e.g. our galaxy. Inserting 
(gH) and ^\ into (|2J leads to the metric: 



A = l- 



2R 



3(l + tj)t 



and the Hubble parameter 



H 



1 1 _ 2 



(46) 



(47) 



where the line-element is given in (|7|). 

The current standard model of the Universe is the flat Friedmann-Lemaitre 
model, which is a Universe model with dust and Lorentz invariant vacuum 
energy; i.e. w = and A > 0. A pedagogical presentation of this model is given 

ininni. 
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For a non zero A there are two cases: A > and A < 0. Integrating (|38|l for 
A > we get 

(48) 



47r cosh [{I + V3A tj - At: 

Taking the series for cosh we get ()44|l in the Umit A — > 0. Eq. H48|l can also be 
written 

P= T-^ (49) 

STTsinh^ Ul+uj)^tj 

Thus, from eq. H21|l the metric ^ is given by 

A = l-^i?2^oth2((l+c.)^t) (50) 
and the Hubble parameter 

((l+c.)^t). (51) 

From 1)40(1 the trajectories are 



/A , .V3A 



( 47r cosh ((1 + cj)\/3At) - 47r\ 3(1+-) 
R = 5 1 (52) 



TTsinlr (1 + w)^ < \ 



= (' '-) (53) 

which are shown for the Friedmann-Lemaitre model, i.e. a; = 0, in figEl Eq. 
I)H5|I then implies that the expansion of the Universe becomes accelerated for 

A > 47r/3, (54) 

or in terms of the Lorentz invariant vacuum energy density: > ^p, '30' . 
Solving H38|l for A < we obtain the solution 

(55) 



247rcos" i^ - ^^^{l+uj)t 
where from H40|l p oc 

This is an oscillating Universe with singularities 
at i? = 0, i = , " — r; n is a half-integer. 

V— 3A(l+(^) ' ° 

The partial differential equation l|25() for M must also give these solutions. 
We now consider this equation for A = 0. Let us set Y = 2m + 2M. In terms 
of Y equation ((25|) for an expanding Universe becomes 

idY , , idY 
i?,_ + (l+^)r._^0, (56) 



11 



which is easily solved by separation of variables. The solution is 

/ 2cJ?J+3c2 

2(m + M) _ l^3(i+^)(ct+ci) 

— 1 7~ — 1 7~ 

R R 
where c is the separation constant and ci and C2 are integration constants. In 
the limit R ^ oo the solution should give eq. (|46|l . this requires that ci — 0. 
Consider now the case where R is finite and t oo. We would in this case, for 
an embedded mass, expect the solution to approach the Schwarzschild solution. 
But here the solution goes towards Minkowski spacetime, indicating that our 
solution represents a pure (m — 0) Universe. Thus, we should set C2 = 0; then 
the separation constant cancels and the solution reduces to (|46|l . This method 
of solving (|25|l cancels out the effect of m. But from this we see that for m = 
equation H25|) gives the required solution. 



3.2 m^O 

We now turn to the case where we have an embedded mass, m 7^ 0. In jl7| 
no discussion was made on solutions of equation H25|l . We consider it for a flat 
Universe with A = 0. 

To find solutions that include the effect of the embedded mass m we shall 
look at approximations. We will find solutions valid close to m and at large 
distances from m. 



IM 1 /M^^ 



We expand the square root 



(2m + 2M)..(2m).(^l + ---— +...), (58) 

thus defining the function F: 

{2m + 2M)^ = (2m) 5 + 2F. (59) 
From ()25|l wc get the following equation for F: 

dF OF OF 

97 + (1 + aR + + ""^^^dR. ^ 

which has two approximations: 

lOF lOF 1 

i?^— + (l+w)(2m)5— «0, 2F<(2m)^, (61) 
Ot or 



1 dt at 1 

i?^— + (l+w)2F— wO, 2F>(2m)^. (62) 

Ot or 

Here H61|l gives an approximate solution for the spacetime close to the embedded 
mass m, and ()62|l gives an approximate solution for the spacetime far from m. 
Both equations are easily solved by separation of variables. 
From equation (|62|l we find the following solution: 

(2m + 2Af)3 = (2m,) ^ + (63) 

3(1 + ujt 
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where we have set the integration constants to zero. The metric is given by 




2R 



3(l + w)i 

This solution is vahd for 2F ^ (2m) ^ which is 

2i?3 



(64) 



9(1+^^)2^2 



> m. (65) 



i.e. where the Shcwarzschild field is much weaker than the background (see l|37(l 
and (|44|l 1 . Taking the hmit i? — > oo this solution approaches the solution H46|) 
obtained for m — Q. Since the solution is only valid far from the embedded 
mass, M is now not a global quantity; i.e. it does not represent the mass of the 
Universe. From H20|l we find the density p of the cosmic fluid with an embedded 



1 (2to)3 
P^'^ ^) = (l + c.)26.t2 + [l + .)A.RH ■ 

The last term gives the deviation from the density H44(l of the cosmic fluid in 
a FRW Universe. We get the geodesies by inserting (|63|l into H34|l . Thus, far 
from m we have 

The solution for = is 

2 

R= (^^V2rn\nt + c^ ti (68) 

where c is the integration constant. We note that for i — > we get R ^ 0. 
Using eq. we see that this solution is valid for: 

i»efS pi = l-_^ (69) 
3(2m)2 



and 



^«e^^ p,^-l--^^. (70) 
3(2m) 2 



For uj the solution is 



Valid for 



R= {^V2^(^ + l]t + ct^ ] . (71) 



^«l rf^l (72) 

,3(2m)i(^) 
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For m = the solutions H68|) and (j71fl reduces to eq. (j45|l with the identification 

Next consider the approximation This equation has solution 

(2m + 2M) ^ = (2m) ^ + ae''"^ , (73) 

where ^ 

a; = r - ^, (74) 

3(l + w)(2m)5 

a is an integration constant and h is the separation constant. This solution is 
valid for 

6x<ln^^ '—. (75) 

a 

The metric function A becomes 

2 



(2m) 



A = l-^ ^ (76) 

For a finite R in the limit i — > oo we get A ^ 1 — 2^, i.e. the solution approaches 
the Schwarzschild solution. Also, for a = 0, which gives M = 0, the second term 
is zero and we are left with the Schwarzschild field. Considering H73|l to first 
order in — this solution reduces to 

m 

M = (2m)5ae^^. (77) 

If we assume M ^ m in l|25|l and ignore the non-linear term completely we can 
write + (2m) = which has solution From ^ and ^ we 

have the geodesies close to m approximated by 

If we write this in terms of x and then insert Z = — = e^^, eq. H78|) can be 
written 

dZ abZ^ - huj{2m)^ Z , , 

— — (79) 

dt {l+u){2m)^ ' 

which is easily integrated to find the trajectories. For dust, a; = 0, the geodesies 



/ 3(2m,)2 3(2m,)2 / ab , , 



are 



2& V (2m) 5 

where Ci is the integration constant. If we demand that R{t = Q) = we get 
Ci = 1. From eq. (|75|l we find that this solution is valid for 

1 / (2m)^ci \ , , 

Thus, for this to be valid in an expanding Universe we must have \/2m ^ a. 
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For Lo ^0 the solution is 



26 \ ac2 + ae 1+" 



Putting R{t=0) — gives the integration constant C2 — ^ — 1 . This solution 

is valid for 

eTf:.^* » (83) 

C2 

By setting a = in (jHOl, or |(H11), we are left with the radial Schwarzschild 
geodesies. 



4 Singular shells in isotropic Universe models 

We are considering a Universe containing, in addition to the cosmic fluid and 
A, energy confined to a surface (or rather an hypersurface) . That is, we are 
dealing with situations where we have a shell of energy where the thickness, d, 
of the shell can be ignored; i.e. mathematically we let the thickness go to zero, 
■)?—!■ 0. This is the thin shell approximation. The energy momentum tensor for 
such a spacetime can be split into three parts. A general Tap can therefore be 
written 

Tap = SapS{y) + T+pe{y) + T-^ei-y), (84) 

here y is an orthogonal coordinate such that ^ = n is a normalized normal 
vector and y = at the hypersurface. 5{y) is a delta function and 9{y) is the step 
function. In the previous section we discussed the embedding of a Schwarzschild 
mass into a cosmological model. Now we wish to obtain the relativistic equation 
of motion for a thin shell in this ambient space-time. From H84() we see that 
the shell contributes with a delta function singularity to the energy momentum 
tensor and thus does not follow geodesies in the background spacetime. In the 
thin shell approximation the energy-momentum tensor Sap of the surface is 
defined as the integral over the thickness of the surface when the thickness •§ 
goes to zero 

± 

Sap = lim / Tapdy. (85) 

•'-2 

4.1 Israel's Formalism: The metric junction method 

To deal with the situation described above we use Israel's formalism In this 
section we give a review. 

The spacetime manifold M is split into two parts, and M^, by a 

hypersurface S. That is, U = A4 with a common boundary: (9A^+ n 
dM" = S. 
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The energy-momentum content of spacetime is coupled to the geometry 
through Einstein's field equations. In the regions and outside the 
hypersurface E we assume that 

G% = 87tT^., (86) 

where + and — means the tensor evaluated in and A4~ , respectively. Thus, 
if T^^ is given, the metrics for the two regions outside S are obtained by solving 
()86|l. The line elements of the two regions are written 

ds^ — g^^dx'^dx'^. (87) 

On the hypersurface S we denote the intrinsic coordinates by and the in- 
trinsic (induced) metric is 

dsl = h,jdCd^\ (88) 

We use Greek indices to run over the coordinates oi {d dimensional) while 
Latin indices run over the intrinsic coordinates of S {d — 1 dimensional). 
The equations for the hypersurface are given by the embeddings (j)"^ : 

^'± = r±{n (89) 

where cc^ are c'" functions with r > 4.* To glue together and Ai^ along 

and g-^ 



the common boundary E we require that the two metrics, 5+ and g , induce 



the same intrinsic metric on E: 

+ dx'^ dx'^ _ dx^ dxt 

We note that the junction is independent of the embeddings x^, which do not 
need to join continuously at the hypersurface. This is an essential property of 
Israel's formalism: We are free to choose coordinates in and A^~, indepen- 
dently. 

Let n be the unit normal to the hypersurface S, which is defined to point 
from M~ to , such that: 



n • n = g^.n'^Ti'' ' ( -1 ^^^^ 

e — +1 gives a spacelike n and thus a timelike hypersurface S. e — —1 gives 
a timelike n and thus a spacelike E. For e = we call the hypersurface a null 
surface^. We shall only consider timelike surfaces, e = +1. 

The induced metric on E can also be written in terms of x^. We have 

h'^ = /i'-' -g^ ■ In terms of the normal vector the induced metric becomes 
= (92) 

*The continuity equations involves the 4th derivative of x^. 

^The formahsm given here breaks down for this case, for a treatment of null surfaces see 

mini EH 
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h'^ defines a projection operator, h^^h" = that picks out the part of a tensor 
that hes in the tangent space of S. 

The extrinsic curvature tensor, K, is essential to the formalism. This tensor 
is a measure of how the hypersurface curves in the surrounding spacetime 
and A4~ . It is defined as the covariant derivative of the normal vector n with 
respect to the connection in and Ai" along the direction of the tangent 
vectors on E. Since n is everywhere normal to S and of constant magnitude, 
the variation of n and thus K will be entirely in the tangent space of S. Also, 
since there is a delta function singularity in the energy momentum tensor at S, 
we get a discontinuity in K at E {S{y) — 0'{y)). The extrinsic curvature tensors 
and in A^"*" and Ai~ , respectively, are thus given by 



± 



(93) 



Note the sign convention. We can express Kij in terms of the Christoffel symbols 



of M"^. Using V(n • ^^e^) = we have 



n, 



(94) 



We see that Kij is symmetric and represents 4-scalars in A^*. In terms of the 



projection operator ^ we have K^^ = -h'^h^Ua./il^ and Kfr = K^, 



In Gaussian normal coordinates the metric is given by 

ds^ = dn^ + hijiC,n)dCd^j, (95) 

where E is located at n = and the induced metric on E is /iy (^*,0) = hij{^^). 
In this case the extrinsic curvature tensor is simply given by Kij = P"^- — 
-i/iy,„ and K'-j = -r\y = -h'^'^kkj^n- 

The Einstein equations on the hypersurface are: 

+2eK^Kjif +e]^g,,{K^ + Ki^K'-^-)f = 8^7^^, (98) 

Integrating H96|l and H97I) across the shell in the thin shell approximation we 
find Sin = Snn = 0. That is, S has no normal components to the shell. We 
have S'^j, = ft." /i^ 5^/3 1^ and 



dx% dxj. ^ |± 
'WW 

Integrating H98|l we arrive at the equation of motion for the surface: 



[K,j]-[K]g,, =8TrS^j. (100) 
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By contraction 

[K..,] = 8^(5,, - \sg,,). (101) 

These equations are called the Lanczos equations and they say that the surface 
energy-momentum tensor is given by the difference in the embeddings of E in 
M.^ . The bracket operation [ ] gives the discontinuity of a tensor at E 

[T]=T+-T-. (102) 

In the same way we define the average { } as 

{r}-i(T++T-). (103) 
We also note two relations between these definitions: 

\TS] = [T]{S} + {T}[S] (104) 

{TS}^{T}{S} + \[T][S]. (105) 

Using [ ] on the Einstein equations (|96(l and (|97|l along with the Lanczos 
equation 1)1 00|) and equation (|l()4|l we find: 

+ [T„„] = 0. (106) 



The discontinuity [T„„] gives the pressure exerted normal to the surface by the 
bulk. Contracting eq. 1)101(1 with Sij we get 

Slra[K'^'] = MS^,S'^ - ]^S^) . (107) 

The contraction Sim.K'^ gives the normal component of the divergence of S 
with respect to the connection on A^='=. 



(108) 



where | denotes covariant derivative with respect to the intrinsic connection on 
E (i.e. the metric connection defined by hij). The normal part is 

n^S''^.^^\^^Si„.K!^. (109) 

Combining (fTUH|l and fTU7|) we can separate SimK^J^ and Si-^K^:^ _13 . In 
a spherically symmetric space-time with the surface consisting of an ideal fluid 
we may separate the time components of and if'™. 

The continuity equation for the surface is 

S- \l + [T^n] = 0. (110) 
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Contracting [Tin] with we get the momentum- flux of the bulk as measured 
by a comoving observer to the surface. Contracting with a spaceUke tangent 
vector, , this term gives the tangential force, in the direction, exerted on 
the surface by the bulk. 

Taking the covariant derivative of the 4-velocity u± in the direction of u-t 
with respect to the connection on we obtain the 4- acceleration a±. As 
viewed from M"^ 



a"e„ |±= Vsw|±= u^^^u^'ei + Ki^u'u^n\^. (Ill) 

The normal component of the 4-acceleration of the shell as viewed from M.^ is 
thus 

a"n„|±=if;t"'«". (112) 

which is in general non-zero. 

Applying { } on H96|l and 19 7|) gives the two constraints: 

- {K^} + {X,„J{if'"} = -IQ^^SimS'^ + S^) - 167r{r„„}(113) 
{K^\i}~{K,,}^~ST,{T,n}. (114) 

4.2 Spherically symmetric thin shells 

We will look at a spherically symmetric shell consisting of an ideal fluid. The 
spherical symmetry implies that we can write the line-element, using proper 
time on the shell, as: 

dsl = -dr^ + Rirfdn'^ (115) 

where R{t) is the expansion factor for the hypersurface, and the proper area is 
given by 

As 47ri?2. (116) 
The energy-momentum tensor for the shell consisting of an ideal fluid is 

Sij = {a+p)uiUj + p gij, (117) 

where a is the mass (energy) density of the surface and p is the tangential 
pressure of the surface. Since we use proper time the 4-velocity of a comoving 
observer is = (1,0,0). The mass, /i, of the surface is 

p = AirR^a. (118) 

The angular coordinates define tangent vectors to the surface. Thus the ra- 
dial coordinates must join continuously on the hypersurface. Hence, the equa- 
tion of the surface is given by 

R±=R{t), (119) 

henceforth omitting the subscripts on R. The line element in the spacetimes 
and when including the effect of a Schwarzschild mass on the background 
is given in ([7||. That is, 
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fi4 = i^-A±dtl-2{l~A±)'!dt±dR + dR^j+R'^dn'^ 

= -dtl+£^^(dR- {1- A±)idt±Y + R^dn^, (120) 



where in the most general case: 

2(m+ + M+) A+ o , , 

A+ = l- ^ +^ (121) 

2(m_+Af_) A_ , , , 

A_ = l--^ ^ i?^ 122) 

R 3 ' ^ ^ 

m_ is the Schwarzschild mass at the center, and to+ = m_ + where is 
the Schwarzschild mass of the shell; A± gives the vacuum energy density in , 
respectively; M± are the mass functions for the Universe in Al*, respectively, 
defined in equation (|19|l . The parameter £± is given in eq. (|32(l . 
The embedding is 

= (t±(r),i?(T),0,0). (123) 
Thus, the 4- velocity, u±, for a comoving observer is 

— ± = -±=^/^ = (4,i?(T),0,0), (124) 
OT dr 

where the dot denotes derivative with respect to the proper time r on the shell. 
In terms of the intrinsic coordinates we have 



86" 



i.e. = (1,0,0). 

From n^Tial^ = 1 and u"na\^ — using = —1 we find the following 

covariant components for the normal vector: 



n 



C±£±\-R.t±,0,0), (126) 
where C. — ±1- 

From the hrr component of the metric junction (|90|l . (or from the line- 
element: dsj. = —dr'^, T is proper time on the shell) we have: 

-A±il ^ 2{1 ~ A±)Hr + £1+ if = (127) 
which is a quadratic equation in t±. Solving this gives: 



(£1- A±)2R^£±JA±+R^ 
t± = 7 ^- ■ (128) 

We make the sign choice by requiring that r and t± are pointing in the same 
direction: i± > 0. Eq. H128|l can be rearranged to give 

t± = ^ , (129) 



il-A±)2R + £±JA±+R^ 
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We sec that this is weU behaved at the horizons A = and crossing the horizon 
such that A < we cannot have stationary sheUs, R ^ 0. Also, we note here 
that t± cannot in general join continuously at the surface when embedding 
in Eq. (|127|l also gives 

^^{£l^A±)^±£±{l-i^^)^ =H±R. + v±. (130) 
at± 

From H33() we see that the first term gives the expansion of the Universe. Thus, 
the second term represents the velocity, v± , of the shell relative to the expansion 
in and , respectively. Then t± is given by the Lorentz factor 

(131) 



1 _ I!± 

^ £± 

From eq. 194(1 we find the angular component of the extrinsic curvature 
tensor, 

Kee = -Rn" = C±R\/a±+R^, (132) 

where we have used (|126f) and H128|l . Thus, from Lanczos eq. H1U1|) we get the 
equation of motion: 



C+^JA++R^-C-\JA^+R^ = ^^A7^Ra. (133) 

This equation generalizes the previous junctions to include the junction for a 
general Universe with an embedded Schwarzschild mass. 

In terms of the Gaussian normal coordinates we have Kgg — —R^\ i.e. the 

(7(7 5 

sign Kgg determines whether the radius of the surface is increasing or decreasing 
in the normal direction. For static spacetimes the sign of Keg is given by C and 
determines the spatial topology IHI IT^ I32 | . If either of the spacetimes are not 
static the sign of Kgg does not in general give the topology; for the classification 
of the junction of two FRW spacetimes see |14| . 
Squaring H133|) we find 



C±\lA^+m = ^^[A+-A^)±^ (134) 

from which one can determine the sign of Kgg. Squaring (|134|l we get the energy 
equation 



The continuity equation H11U|) gives for a comoving observer: 

a^-2{a+p)'^ + [T^,,n'^u''], (136) 

or in terms of /i we have 

fi = -SnpRR + 47ri?2 [T^^n'^u''] . (137) 
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The first term is due to the tangential pressure of the surface. The second 
term, the area of the shell times the discontinuity in the 4-momentum, can be 
interpreted as the mass gathered on the shell from the surroundings due to the 
motion of the shell relative to the cosmic fluid. The force needed to move on 
the surface is given by 

{a+p)u,\iu' = -{5\ + u,u')p^i - u.iT^^n'^u''] - [T^.n^^], (138) 

which for a comoving observer is zero, i.e. following the intrinsic geodesies of 
the surface. 

The energy-momentum tensor is given in equation l|15|l . where (eq. ©) 
is the 4-velocity of the cosmic fluid. Contracting with the normal and the 

4-velocity of the surface gives: V^n^l^ = —9± ^(1 — A±)ii± — Rj, and 

V^u^l^ = —i±. Thus, when including the vacuum energy, the components of 
the energy-momentum tensor at the surface becomes: 



w^n* = (p±+P^/)(i?-(i-A±)^4) +pff-^ (139) 

T^.n^u-'l^ = -C±{p±+p%')i±(^R-{l-A±)n±y (140) 

where Pcf is the pressure of the cosmic fluid and p the energy density. For a 
shell moving with the expansion of the Universe we have R = (f^ — A±)i and 
i± = 1. This gives T^.^nt^n" |±= p^^ - ^ and Tf.^n'^u'' |±= 0. 

For an equation of state pcf — ujp, M satisfies eq. H25(l . Thus, the derivative 
of M with respect to the proper time r on the surface may be written: 



M± = ^(4-(l + c,±)(4-A±)-^7?) (141) 

M± = ^(i?-(l+^±)(4-^±)^4). (142) 

We will in the following restrict our discussion to = 0. We see that if the shell 
is moving with the expansion this is zero, M = 0. Also, we find that Af > 
for a shell expanding faster than the Universe and M < if the shell is moving 
slower than the expansion of the Universe. Using eq. (|142|l with = the 
normal and mixed components of the energy momentum tensor (|139|l and H14()(l 
can be written: 



M±t± 



T,.n^u^\- = ~C±^. (144) 

Furthermore, M can be written in terms of v± and p± by using equations H20|l . 
(fnn|l and jnill: 

M^^'-^^££±. (145) 
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This gives the rate of change of the mass /i, 



/i = -SnpRR - [ivATTR^p] (146) 

where we have set C± = ^1- For a surface with vanishing tangential pressure, 
p = 0, this has the same form as eq. 2.13 in |15| . 

Consider now the time component of the field equations for the surface. 
Combining equations HlUt)|) and (|107|l in the case where the shell consist of an 
ideal fluid (|117|l we find in general 

where the superscripts ± correspond to the ± in the equation, respectively, [a is 
the energy density of the surface) . For a comoving observer K^^ gives the proper 
(normal) acceleration, i.e. K^^ — d^nf^. Hence, the terms on the right in eq. 
(|147|l represent the forces on the shell making it deviate from geodesic motion 
in tW"*" and A^~, respectively. The first term is the difference in the normal 
force exerted on the shell from Al*, the second term is due to the tangential 
pressure on the surface while the last represents the self gravity of the shell. 

We shall consider eq. 1)147(1 for a flat Universe with dust and vacuum energy, 
i.e. Pcf = and £ = 1; also we put K± = A. Differentiating eq. 1134(1 with 
respect to t, then using ifT^ . ifT^ . ifT^ . lfTT7)l . (fH^ and ifT^ we find 



Ml 

dM± 
OR 



± 2TTa ± Snp, 



(148) 



where 



R 



m± 

'W 



M± A± p 



For a test shell, Krr = 0, we recover the Friedmann eq. ((35(1 . 



(149) 



5 Conclusion 

In this paper we have discussed the embedding of a Schwarzschild mass into 
a cosmological model using "curvature" coordinates. Extending the work by 
Gautreau jl7[ we have found approximate solutions to eq. I(25() giving the mass 
function M for the Universe explicitly, and we have solved for the radial geodesies 
outside the embedded mass. We have considered spacetime close to m, where our 
solutions go towards the Schwarzschild spacetime. And far from m, our solutions 
approach the FLRW Universe models. In particular we have presented solutions 
for a flat Universe with vanishing cosmological constant and an equation of state 
p — up, Lu — constant. Without the embedded mass we have given solutions 
with and without cosmological constant for the equation of state p = up. 

Using the Gautreau metric we have generalized Israel's formalism to singular 
shells in a Sch/FLRW background. For an arbitrary equation of state for the 
surface the equations governing its motion are ((135(1 , 1(137(1 and ((148(1 . Equations 
((135(1 and p37(l are given for the general case, while we have only considered 
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the acceleration (|148ll for a flat pressurefree Universe model. An interesting 
further development would be to extend this to three-dimensional branes in a 
five-dimensional spacetime. 

An application of the results obtained in this paper is to study the evolution 
of cosmological voids, see TS", TB", and references therein. In the collapse 
of a positive perturbation leading to the Einstein-Straus model is discussed, 
i.e. the expansion of the voids are comoving. In |15l [T^ a less dense region is 
considered. This less dense region will expand faster than the outer region and 
numerical simulations show that a thin shell is formed. It would be interesting 
to investigate if this problem can be analyzed analytically. Also, it would be 
interesting to apply our formalism to slowly rotating shells, generalizing the 
traetments in ^11201 to shells that need not be comoving with the cosmic fluid. 
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